The Large-Scale Polarization of the Microwave Foreground 
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Most of the useful information about inflationary gravitational waves and reionization is on large 
angular scales where Galactic foreground contamination is the worst, so a key challenge is to model, 
quantify and remove polarized foregrounds. We use the Leiden radio surveys to quantify the polar- 
ized synchrotron radiation at large angular scales, which is likely to be the most challenging polarized 
contaminant for the WMAP satellite. We find that the synchrotron E- and B-contributions are equal 
to within 10% from 408 — 820MHz with a hint of i5-domination at higher frequencies. We quantify 
Faraday Rotation & Depolarization effects and show that they cause the synchrotron polarization 
percentage to drop both towards lower frequencies and towards lower multipoles. 
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I. INTRODUCTION 

CMB polarization and its decomposition into E and 
B modes is a topic of growing importance and interest 
in cosmology [1]. In the era of WMAP [2], a key issue 
is to estimate the contribution of Galactic foregrounds 
(more specifically, polarized synchrotron emission) at the 
large angular scales. Unfortunately, these large scales are 
also the ones where polarized foreground contamination 
is likely to be most severe, both because of the red power 
spectra of diffuse Galactic synchrotron and dust emission 
and because they require using a large fraction of the sky, 
including less clean patches. The key challenge in the 
CMB polarization endeavor will therefore be modeling, 
quantifying and removing large-scale polarized Galactic 
foregrounds. 

Unfortunately, we still know basically nothing about 
the polarized contribution of the Galactic synchrotron 
component at CMB frequencies [3-9], since it has only 
been measured at lower frequencies and extrapolation is 
complicated by Faraday Rotation. This is in stark con- 
trast to the CMB itself, where the expected polarized 
power spectra and their dependence on cosmological pa- 
rameters has been computed from first principles to high 
accuracy [10-13]. 

This is the topic of the present proceeding. We will em- 
ploy polarization sensitive radio surveys to further quan- 
tify the polarized synchrotron radiation, which is likely 
to be the most challenging contaminant in the polari- 
zation maps expected from the WMAP satellite [14,15]. 
This proceeding is organized as follows: in section Sec- 
tion II, we review the basics of the synchrotron emission, 
as well the problems involved with extrapolations from 
lower to higher frequencies. We present our results as 
well as discuss our conclusions in Section III. For more 
details about this analysis consulte [16] 



II. OUR KNOWLEDGE OF SYNCHROTRON 
EMISSION 

The Galactic InterStellar Medium (ISM) is a highly 
complex medium with many different constituents inter- 
acting through a multitude of physical processes. Free 
electrons spiraling around the Galactic magnetic field 
lines emit synchrotron radiation [17], which can be up 
to 70% linearly polarized (see [18,19] for a review). 

The power spectrum Ct of synchrotron radiation is 
normally modeled as a power law in both multipole £ 
and frequency is, which we will parametrize as 



STf(u) = A[-) with A 



P+2 



(1) 



where 6T e = [£(£ + 1)Ci/2tt} 1 / 2 . This definition implies 
that Cf oc 1$ for I 3> 1 and that the fluctuation ampli- 
tude oc v a . The standard assumption is that the unpo- 
larized intensity has a ~ —2.8 with variations of order 
0.15 across the sky [20] - see also [21-23]. 

As to the power spectrum slope (3, the 408 MHz 
Haslam map [24,25] suggests (3 of order -2.5 to -3.0 down 
to its resolution limit of ~ 1° [26-29]. A similar analysis 
done on the 2.3 GHz Rhodes map of resolution 20' [22] 
gives 13 = -2.92 ± 0.07 [30] (flattening to /3 « -2.4 at 
low Galactic latitudes [9]). 

For the polarized synchrotron component, our ob- 
servational knowledge is, unfortunately, not as com- 
plete. To date, there are measurements of the polar- 
ized synchrotron power spectrum obtained basically from 
three different surveys [31]: the Leiden surveys [32,33], 
the Parkes 2.4 GHz Survey of the Southern Galactic 
Plane [34,35], and the Medium Galactic Latitude Sur- 
vey [36-38]. These measurements exhibit a much bluer 
power spectrum in polarization than in intensity, with (3 
in the range from 1.4 to 1.8 [3-9]. These results are usu- 
ally taken with a grain of salt when it comes to their im- 
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plications for CMB foreground contamination, for three 
reasons: 

1. Extrapolations are done from low to high latitudes; 

2. Extrapolations are done from low to high frequencies; 

3. Much of the available data is undersampled. 

The Leiden surveys extend to high Galactic latitudes 
and up to 1.4 GHz but arc unfortunately undersampled, 
while the Parkes and the Medium Galactic Latitude Sur- 
veys only probe regions around the Galactic plane - see 
[16] for more details. In the following three subsections, 
we will discuss these three problems in turn. 

A. The Latitude Extrapolation Problem 

There is a well-know empirical result that shows that 
whereas the unpolarized synchrotron emission (at MHz 
range) depends strongly on the Galactic latitude, the 
polarized component is approximately independent of 
Galactic latitude (see, e.g., [35]). The usual interpreta- 
tion for this very weak latitude dependence of polarized 
synchrotron radiation is that the signal is dominated by 
sources that are nearby compared to the scale height 
of the Galactic disk, with more distant sources being 
washed out by Depolarization (to which we return in the 
next subsection). As a result, having well-sampled polar- 
ized maps off the galactic plane at the same frequencies 
would not be expected to affect our results much, since 
they would be similar to those in the plane. This issue, 
however, deserves more work as far as extrapolation to 
CMB frequencies is concerned: the latitude dependence 
may well return at higher frequencies as Depolarization 
becomes less important, thereby revealing structure from 
more distant parts of the Galactic plane. In this case, ex- 
trapolating from an observing region around the Galactic 
plane to higher latitudes may well result in less small- 
scale power in the angular distribution. 

B. Frequency Extrapolation Problem 

It is important to point out, that Faraday Rotation 
(see, e.g., [39]) can only change the polarization angle 
and not the polarized intensity P {P=\/Q 2 + U 2 ). The 
fact that we do see structure in P that is not correlated 
with a counterpart in intensity T implies that part of the 
radiation has been depolarized [40]. Depending on the 
frequency and beamwidth used, Depolarization can play 
an important role in polarization studies of the Galac- 
tic radio emission [33] - for more details, see Cortiglioni 
and Spoelstra [41]. Because of the complicated interplay 
of these mechanisms, we should expect both the ampli- 
tude and the shape of the polarized synchrotron power 
spectrum to change with frequency. 



Table 1 Normalization &z Spectral Index( a ) 



v A E 


Pe 


A B 


Pb 


(GHz) [mK 2 ] 




[mK 2 ] 




0.408 5.5 


-0.5 


5.7 


-0.4 


0.465 5.4 


-1.0 


5.4 


-0.5 


0.610 5.1 


-1.0 


5.1 


-0.8 


0.820 4.5 


-1.5 


4.6 


-1.8 


1.411 3.9 


-1.9 


3.6 


-2.6 



(a) All fits are normalized at £=50, i.e., 8Tj = A(£/5O) 0+2 . 



C. Incomplete Sky Coverage and the Undersampling 
Problem 

For the case of undersampling in the Leiden surveys, 
some authors have overcome this problem by doing their 
Fourier analysis over selected patches in the sky where 
they believe the average grid space in the patch is close 
to the map's beam size, so that they can apply a Gaus- 
sian smoothing on it - this is well explained and illus- 
trated in [7]. Fortunately, we can eliminate this problem 
by measuring the power spectra with the matrix-based 
quadratic estimator technique that has recently been de- 
veloped for analyzing CMB maps [42-44]. Although the 
undersampling and partial sky coverage results in un- 
avoidable mixing between different angular scales I and 
polarization types (E and B), this mixing (a.k.a. leak- 
age) is fully quantified by the window functions that our 
method computes [43] and can therefore be included in 
the statistical analysis without approximations. Specifi- 
cally, we compute the six power spectra (Cf = T, Cf = 
E, Cf = B, Cj E = X, Cj B = Y, Cf B = Z) so that the 
much discussed [43,12,45-49] E — B leakage is minimal 
[44]. 

III. RESULTS & CONCLUSIONS 

We employed only the Leiden surveys [32,33] for our 
analysis. The observations done by Brouw and Spoelstra 
covered almost 40% of the sky extending to high Galac- 
tic latitudes. Using the same instrument, they observed 
the polarized Galaxy in Q and U in five frequencies from 
408 MHz up to 1.4 GHz and with angular resolutions 
from 2.3° at 408 MHz up to 0.6° at 1.4GHz. Unfortu- 
nately this data was also undersampled, making it diffi- 
cult to draw inferences about its polarized power spec- 
trum. 

Using matrix-based quadratic estimator methods 
[43,44], we measure the power spectra from the Leiden 
surveys, obtaining the following key results: 

1. Our analysis was performed using 10 multipolc 
bands of width A£ = 10 for each of the six polariza- 
tion types (T, E, B, X, Y, Z), thereby going out to 



2 



(, = 100. We used the Haslam map for the unpolar- 
ized component T, scaled and smoothed to match 
Leiden's five different frequencies, and assuming a 
|6| = 25° Galactic cut. The best fit normalizations 
A and slopes j3 for E and B are shown in Table 1. 
The values of (3 are consistent with previous anal- 
yses [3-9], showing that the slopes get redder as 
frequency increases. For all Leiden surveys, the 
X and Y power spectra are found to be consistent 
with zero - the 2.4 GHz Parkes survey had a similar 
finding for X [9]. These are not surprising results: 
if Faraday Rotation makes the polarized and unpo- 
larized components to be uncorrelated, it is natural 
to expect that X, Y—0. However, at the CMB fre- 
quencies (where the effects of Faraday Rotation & 
Depolarization are unimportant) this should not be 
the case. 

2. To study the frequency dependence, we average the 
10 multipole bands of the Leiden power spectrum 
measurements together into a single band for each 
polarization type to reduce noise. From these re- 
sults, we fit the average frequency dependence (for 
the 25° cut data) as a power law as in equation (1) 
with slope cxe = —1.3 and as = —1.5 for E— and 
B— polarization, respectively. 

3. An interesting question about polarized fore- 
grounds is how their fluctuations separate into E 
and B. Although many authors initially assumed 
that foregrounds would naturally produce equal 
amounts of E and B, Zaldarriaga [46] showed that 
this need not be the case. Early studies [5,9] have 
indicated that E w B at 2.4 GHz in the Galactic 
plane. However, these analyses used Fourier trans- 
forms and spin- 2 angular harmonic expansions, re- 
spectively, without explicitly computing the win- 
dow functions quantifying the leakage between E 
and B. We therefore perform a likelihood anal- 
ysis of the Leiden surveys specifically focusing on 
this question, and including an exact treatment of 
the leakage. The likelihood analysis of the data is 
done with two free parameters corresponding to the 
overall normalization of the E and B power spec- 
tra, and assuming that they both have the same 
power law shape given by the slopes 0e from Ta- 
ble 1. We obtain that the synchrotron E- and B- 
contributions are equal to within 10% from 408 to 
820 MHz, with a hint of _E-domination at higher 
frequencies. One interpretation is that E > B at 
CMB frequencies but that Faraday Rotation mixes 
the two at low frequencies. 

4. Faraday Rotation & Depolarization effects depend 
not only on frequency but also on angular scale - 
they are important at low frequencies (y < 10 GHz) 
and on large angular scales. Therefore, we must 
take into account Faraday Rotation & Depolariza- 
tion effects when extrapolating radio survey results 



from low to high galactic latitudes and from low to 
high frequencies. 

5. We detect no significant synchrotron X cross cor- 
relation, but Faraday Rotation could have hidden a 
substantial correlation detectable at CMB frequen- 
cies. 

6. Combining the POLAR [50,51,16] and radio fre- 
quency results, and the fact that the i?-polarization 
of the abundant Haslam signal in the POLAR re- 
gion is not detected at 30 GHz, suggests that the 
synchrotron polarization percentage p at CMB fre- 
quencies is rather low (p <20%). 

Experiments such as polarized WMAP and Planck will 
shed significant new light on synchrotron polarization 
and allow better quantification of its impact both on 
these experiments and on ground-based CMB observa- 
tions. 
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